
control = readMesh('cube.obj');
deform = readMesh('shear3x3.obj');
mesh = readMesh('cube3x3.obj');

rc = control.v * 0.2 + 0.4;

rx = rc(:, 1);
ry = rc(:, 2);
rz = rc(:, 3);

echo = 'mesh loaded'

A = assemblyMatrix(mesh);

echo = 'matrix assemblied'

[o, b] = assemblyBoundaryVector(mesh, deform); 

echo = 'boundary assemblied'

theta = A \ b;

echo = 'matrix solved'

u = evaluateInteriorPoint(rx, ry, rz, mesh, theta);

echo = 'interior points evaluated'

plotDeformation(mesh, rx, ry, rz, o, u);

